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Abstract. We show that well-separated quantum superposition states, measure- 
ments of strongly nonlinear observables, and quantum dynamics driven by anomalous 
diffusion can all be achieved for single atoms or molecules by imaging spontaneous 
photons that they emit via resonance florescence. To generate anomalous diffusion 
we introduce continuous measurements driven by Levy processes, and prove a number 
of results regarding their properties. In particular we present strong evidence that 
the only stable Levy density that can realize a strictly continuous measurement is the 
Gaussian. 
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1. Introduction 

The measurement of the position of a particle is perhaps the most basic quantum- 
mechanical measurement, but it poses many subtleties. The archetype for a position 
measurement is the Heisenberg microscope — the imaged detection of light scattered from 
the particle. Here we consider photodetection of the light emitted by a single two-level 
atom. If the photons first pass through a lens designed to form an image of the atom 
before detection, then the photodetections provide information about the location of 
the atom. The detections therefore modify the spatial wave function of the atom. This 
situation was first considered by Holland et al. [1], as a tool for the efficient simulation 
of atomic decoherence due to spontaneous emission. As applications, they considered 
spatial flights correlated with quantum jumps, and Levy statistics due to trapping in 
potential wells of an optical lattice [2]. 

Under the proper conditions, this type of imaged photodetection can lead to the 
standard form of a continuous position measurement [3], which gives rise to Gaussian 
projection noise [4]. However, as we will discuss in detail, it is the diffraction pattern 
of the imaging system that gives the form of the collapse in the Heisenberg microscope. 
Imaging systems often involve hard-edged apertures and thus long-tailed (non- Gaussian) 
diffraction patterns. Such long-tailed distributions lead in the context of random walks 
to exotic Levy noise — so-called "anomalous diffusion" [5, 6, 7]. Another situation in 
which one can naturally tailor the form of the collapse would be the detection of 
resonance fluorescence in the presence of a magnetic field gradient [8]. In this case it is 
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the Lorentzian line shape that enters as the collapse function. Thus, in moving towards 
experiments to realize continuous position measurements via imaging, it is important to 
understand the impact of long-tailed collapse operators, and under what conditions there 
may be qualitative modifications to the quantum noise and the conditioned dynamics 
of the continuously observed atom. 

Our purpose here is threefold. In Sec. 2 we examine the effect of imaging on 
the wave-function of a single atom, and show that this can be used to prepare non- 
classical states, and realize highly nonlinear measurements. In Sec. 3, in preparation 
for exploring how imaging can be used to generate anomalous diffusion, we examine 
continuous measurements that contain Levy noise. We are able to show that there are 
no truly continuous measurements that are driven only by stable Levy processes [9, 10], 
with the sole exception being the Gaussian. It turns out that the underlying reason 
for this is a result that we prove in Sec. 4: repeated measurements whose errors are 
given by the stable Levy distributions will eventually collapse the wave-function to a 
Gaussian, just like the usual Gaussian measurements. In Sec. 3 we also discuss how 
truly continuous measurements can contain Levy noise, by combining a stable Levy 
process with a Poisson or Gaussian process. Finally, in Sec. 4 we use the above results 
to examine how imaging can be designed to make measurements that induce anomalous 
diffusion, an effect that, in this case, is quantum-mechanical in origin. 

2. Position measurements via photodetection 

We consider the motion of a single two-level atom (or molecule) along a single direction, 
which we will refer to as the z-axis. We illuminate the atom with a resonant laser 
traveling along the x-axis, and place two lenses or mirrors on opposite sides of the 
atom that image the emitted photons in the xz-plane. This configuration is depicted in 
Fig. 1. In addition to forming an image of the atom at the detector, the mirrors have 
a mask over them, which can apply an angle-dependent phase, as well as having an 
angle- dependent absorption profile. This mask introduces an aperture for the imaging 
optics. 

To describe the angular dependence of the aperture, we use spherical polar 
coordinates 9 and <fi. As usual 9 is the angle to the z-axis, and <fi is the angle in 
the x?/-plane. The center of each lens is on the y-axis, and thus given by 9 = n/2. We 
parameterize the distance from the center of each lens by £ = cos 9 G [—1,1]. We will 
denote the aperture transmission function of the aperture, 0), and for simplicity we 
will take this to be independent of in a fixed window A<p (centered at the center of the 
lens), and to be zero otherwise. We note that since the aperture may include a phase 
mask, t can be complex. 

The basic analysis of the above imaging setup has been performed in Refs. [1, 3], 
and the relevant results are as follows. Upon detecting a single photon at location z = a 
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Figure 1. Spontaneous emissions from a single atom or molecule are imaged by 
mirrors. Each mirror has a mask whose transmission profile can be modified. 



on the image, the wavefunction of the atom, ip(z), undergoes the transformation 

where M is the required normalization, and the "collapse" operator A(z) is 

A(z) = Jm^dt (2) 
Here k is the wave number of the emitted photon, and 



x(o=mvW), (3) 

where = (3/4) sin 2 (9) is the dipole angular emission function for the photon. 

By choosing the aperture transmission appropriately, x can be any function satisfying 
|x(OI 2 — /(0) with the restriction that x is zero outside the aperture capture region. 
The fraction i] of photons captured by the optics is given by the integral of \x\ 2 over 
the capture region. An emitted photon that is not detected transforms the atomic wave 
function via the collapse operator B(z — a), where B(z) is given by Eq. (2), with x 
replaced with x' = V7(l ~~ t). Since the photon is not detected, the density matrix of 
the system is given by averaging over all the possible image locations a. 

The power of the above imaging setup comes partly from the ability to engineer the 
collapse (or measurement) operator A(z) by selecting an appropriate aperture profile, 
Since is zero outside [—1, 1], for a given A(z) we can find the required profile 
by using the inverse Fourier transform: 

t(0 oc -== / A(z)e- ikz ^dz. (4) 

The other key component of imaging is the ability to engineer the mixing of optical 
modes prior to detection. We now present three examples that illustrate how these two 
components can be exploited. 
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Figure 2. Aperture amplitude transmission profiles for preparing (a) a spatially 
separated superposition, and (b) a state with an infinite position variance. 



2.1. Application: Preparing spatially separated superpositions 

The imaging configuration described above can be used to prepare an atom in a 
superposition of two spatially separated states. This is achieved by chosing the aperture 
of the imaging system so that the collapse function A(z) is the sum of two well-separated 
Gaussians. Let a 2 be the variance of both Gaussians, and L 3> a their separation. An 
aperture function that does this is shown in Fig. 2. This aperture covers about 29° 
of arc in 9 (i.e. £ G [—1/4,1/4]), and gives well-separated Gaussians with a = 1.5A 
and L = 15A. Here A is the wavelength of the emitted photon. The separation is 
easily increased by increasing the number of oscillations in the profile. The fraction of 
photons captured by the aperture, assuming the two- mirror configuration in Fig. 1, and 
A<p = 29°, is approximately 1/186. The aperture size we have chosen is conservative - 
using two mirrors with 120° of arc as the focussing elements, capture rates greater than 
1 in 20 should be feasible with this aperture profile. 

One now prepares the atom in a broad Gaussian state, centered at z — 0, by 
cooling it to the ground state of a harmonic trap using standard techniques [11]. One 
then illuminates the atom with a laser pulse to prepare it in the excited state, and 
allows it to emit a photon. This preparation process is repeated until the emitted 
photon is detected. (One must repeat the preparation because a single undetected 
emission destroys the spatial coherence.) This prepares the atomic wave function in 
a superposition of two spatially separated wave-packets, with separation L and widths 
~ a. The relative heights of the two localized wave-packets is determined by the location 
of the detected photon. The closer the photon to z — 0, the more equal the weighting 
of the two packets. Naturally one requires that the probability of the smaller packet, 
P s , should not be too small. If the ground state of the trap is broad compared to L, 
a numerical calculation shows that when a photon is detected, the probability that it 
gives P s > 1/3 is approximately 95%. Furthermore, the superposition can be made 
arbitrarily symmetric with a sufficient number of repetitions of the preparation process. 
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2.2. Application: Measuring nonlinear observables 

Measurements of nonlinear observables have applications in state-preparation [12] and in 
detecting signatures of quantum motion, among others. The usual method of measuring 
nonlinear observables (e.g. nonlinear functions of position, z), is to couple a probe to the 
system via a nonlinear interaction proportional to the observable. Generating effective 
nonlinear interactions is not easy, however. Feasible methods exist to generate effective 
interactions, but generating higher powers becomes increasingly more difficult [13]. 
We show now that imaging can be used to engineer a measurement of the highly 
nonlinear function, \z\, without the need for a nonlinear interaction. This measurement 
can be used to generate spatially separated superposition states, in a similar manner 
to that of z 2 measurements [12], with the advantage that it is linear everywhere but at 
z = 0. This linearity means that the wave-packets in the resulting superposition will be 
very nearly Gaussian squeezed states. 

To realize a measurement of \z\ we make the aperture fully transmitting. This 
allows one, in principle, to capture the majority of emitted photons. This time we use 
the optics to first collimate the light. We then split the light into two beams with a 
50-50 beam splitter, invert one of them about z = 0, and recombine them. Finally we 
focus the beam to an image, which is now a superposition of the original image and 
its reflection about z — 0. Since we cannot now tell z from —z, the detected photon 
only contains information about \z\. This measurement can be made quasi-continuous 
by repeatedly exciting the atom and detecting the emitted photons. The method can 
also be extended to more complex measurements, erasing information about multiple 
intervals of the values of z. 

3. Anomalous diffusion and Levy measurements 

All continuous measurements that have been realized to date on quantum systems 
generate Gaussian or Poissonian noise. That is, in each infinitesimal time step dt, 
the error in the measurement record, and the noise induced in the system, is either the 
Gaussian (Wiener) increment dW, with variance dt, or the discrete jump of a Poisson 
process. These noise processes are ubiquitous in nature because of the central limit 
theorem. There are, however, a whole range of more exotic noise processes that do 
not have these statistics, and find applications in wide range of phenomena [5, 6]. A 
remarkable class of these are the symmetric stable Levy processes, L a (t), indexed by 
the continuous parameter a £ (0,2] [9, 10]. We will denote the random infinitesimal 
increment of one of these processes by dL a . When a = 2, the Levy process is just the 
Wiener process. All other stable Levy processes break the central limit theorem because 
they have infinite variances. 

The "stable" Levy processes are so-called because the sum of two or more of their 
increments has the same probability density as each individual increment (except that 
the width of this density is larger by a factor of 2 1 / a ). This is familiar for Gaussian 
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random variables (the sum of two Gaussian deviates is also Gaussian), but is usually 
not true for other probability densities. In fact, if a probability density satisfies the 
conditions for the central limit theorem, and is not Gaussian, then the sum of multiple 
samples from this density must differ from the original density, since ultimately this 
sum must have a Gaussian density in the limit of many samples. 

When two increments of the (Gaussian) Wiener process are added together, the 
width (uncertainty) of the sum is v2 times that of each increment. As stated above, 
for the stable Levy process L a this factor becomes 2 1 /". This factor determines how 
the width (the uncertainty) of the process scales with time. The result is that the 
process L a (t) has a width proportional to t l l a . For Gaussian noise this reduces to the 
familiar scaling: the standard deviation increases as y/t. This is the usual behavior of 
diffusion, for example the diffusion of a pollen grain undergoing Brownian motion. If the 
uncertainty of a dynamical system scales instead as for a < 2, then this behavior 
is referred to as anomalous diffusion [7]. The case a = 1 is the Cauchy process, whose 
increment, dC = dL±, has the Lorentzian probability density 

p{dc) = W+WT (5) 

The uncertainty of the Cauchy process scales at t. For readers new to Levy processes, 
realizations of the Wiener and Cauchy processes are displayed in Fig. 3. A realization, 
or "sample path" , of the Cauchy process is quite distinct from the Wiener process. An 
introduction to Levy processes may be found in [9, 10]. With this background out of 
the way, we ask how stable Levy processes might be realized by continuous quantum 
measurements. 

3.1. Absence of continuous measurements with stable Levy noise 

A continuous measurement of position, X, is characterized by the error of the 
measurement result in an infinitesimal time-interval dt. Specifically, if de is the random 
variable describing this error, then the increment of the measurement result obtained in 
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the time-interval dt is given by 



dr = (X)dt + 7 de, 



(6) 



where 7 is a constant that fixes the overall size of the error. Gaussian measurements 
have de = dW [3, 14, 15]. Little is yet known as to what classes of exotic Levy processes 
can appear as the errors in continuous measurements [16]. We now show that continuous 
measurements can never have errors given purely by stable Levy processes, with the 
sole exception being the Gaussian. 

The amount of information that the measurement obtains about X in a short time 
interval At is determined by the width of the probability density of the measurement 
result in each time interval. The measurement result is given by dividing dr by dt, and 
is the mean value of the measured quantity plus the random error: 



The amount of information extracted is determined by the width of the probability 
density of de/dt; if the error has a large "variance" (width) the observer gains little 
information, and vice versa. In particular, the observers probability density for X (the 
diagonal elements of the density matrix in the basis of X) after the measurement is 
given by multiplying her initial probability density by the probability density for de/dt, 
shifted by the measurement result, dr/dt [3]. 

For a sequence of measurements to produce a valid continuous measurement in 
the continuum limit, the amount of information extracted (the change to the observers 
state-of-knowledge) must scale in the right way as we reduce the time interval between 
the measurements, At. In particular, the average change to the observer's state-of- 
knowledge in a fixed time interval, T, must remain the same as we reduce the time-step 
and thus increase the number of measurements in the interval. For this reason we can 
show immediately that Eq. (6) cannot describe a continuous measurement if de = dL a 
for a < 1. Consider first the case when a — 1. Because the width of dL\ scales as dt, we 
see that the width of de/dt = dL\/dt is independent of dt. Since it does not decrease with 
dt, as we increase the number of measurements, and take the limit dt — > 0, the amount 
of information extracted by the observer tends to infinity, resulting in an instantaneous 
collapse of the wave-function. This will also be true for a < 1. 

Showing that continuous measurements cannot be realized for 1 < a < 2 is 
not so simple, since in this case the width of the probability density for the error, 
de/dt = dL a /dt, diverges as dt — > 0, and thus the extracted information tends to zero 
as required. To determine whether a continuum limit exists for these measurements we 
begin by way of a numerical example before we proceed with an analytic calculation. 
We start with a state of knowledge for X given by a Gaussian with unit variance, and 
simulate a sequence of measurements for a = 1.5. This involves repeatedly multiplying 
by the probability density for L a , which must be obtained from the Fourier transform 
of its characteristic function [9]. We simulate four sequences of measurements, in which 
each sequence has half the time-step of the one before. We then examine how the change 



dr 
~dt 



W+7 



de 
dt' 



(7) 
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Figure 4. The scaling of the average reduction in an observer's uncertainty, for a 
sequence of stable Levy measurements, as the time step is decreased, and thus as the 
number of repetitions is increased. The horizontal axis gives the number of repetitions. 
The straight line gives the best (logarithmic) fit. 

in the width of the observer's state-of-knowledge scales with the time-step. We obtain 
this change by averaging each sequence over four thousand realizations. The results for 
the four sequences are displayed on a log-log plot in Fig. 4. From this we see that the 
reduction in the observers uncertainty increases as a power of the inverse time-step (the 
exponent, or slope in the plot here, is consistent with 1/3), and thus will reduce the 
observer's uncertainty to zero in the limit dt — > 0. A sensible continuum limit would 
require a flat slope (flat scaling exponent) in such a plot as this, as happens in the 
Gaussian case a = 2. Again, this means that the amount of information extracted per 
unit time becomes asymptotically constant as dt — > 0. 

To test this for all a G (1,2), we can compute the scaling exponent as a function 
of a. To perform this calculation numerically, we can use a deterministic method 
that is more efficient than the Monte Carlo method of Fig. 4. The idea is to take a 
narrow, Gaussian initial state, and compute the uncertainty reduction that results from 
applying the measurement operator centered about a particular measurement result. We 
then average the reduction in uncertainty over all measurement results by performing a 
numerical integral. 

The above numerical procedure for calculating the scaling exponent as a function 
of a is rather cumbersome, since there is not a general expression for the measurement 
operator for arbitrary a. Remarkably we can bypass this procedure, and obtain a simple 
analytic expression for this scaling exponent, using the following insight: if the initial 
state is Gaussian, the application of the measurement function (Levy density) leaves 
the final state in a Gaussian distribution, and the reduction of the variance behaves, on 
average, as if the measurement function were itself Gaussian. This is not an obvious 
statement, and we will justify it in Sections 4.1, 4.2, and Appendix A. This result 
means that the reduction in the variance due to the collapse is simple for every Levy 
measurement: it is proportional to the inverse square of the width of the distribution 
for the error). Now that we know how the reduction in the variance depends on 
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Figure 5. The dependence of the power-law scaling exponent (i.e., the slope of the 
line in Fig. 4) as a function of a. The solid curve is the analytic result (2/a — 1), 
and the points indicate values obtained by a numerical calculation of the uncertainty 
reduction due to projection with the corresponding Levy distribution. 

the width, and we know how the width of the Levy distribution scales with dt, we 
can easily determine the scaling of the reduction in the variance. Specifically, if the 
measurement error is given by dL a /dt, then the square width of this distribution is 
w 2 (dt) = dt 2 l a ~ 2 . So the reduction in the variance for a single measurement for duration 
dt is AV(dt) oc —l/(w 2 ) = —dt 2 ~ 2 l a . If we now make N measurements, each with 
duration dt/N, then the total reduction is 

AVtot = NAV(dt/N) oc -N(dt/N) 2 ~ 2/a = N 2/a ^AV(dt). (8) 

The scaling exponent we seek is therefore 2/a — 1. 

Fig. 5 shows the dependence of the power-law scaling exponent (i.e., the slope 
of the line in Fig. 4) as a function of a. The analytic result 2/ a — 1 is shown 
along with numerical tests at the specific points indicated. Again, only a scaling 
exponent of zero can correspond to a sensible continuous limit of a sequence of Levy 
measurements. Thus we see that, even in the range 1 < a < 2, the uncertainty 
decreases too quickly to correspond to a continuous measurement, as the continuous 
limit will lead to instantaneous collapse of the state — essentially, corresponding to a 
divergent measurement strength. Though we have only analyzed a Gaussian state, 
a valid continuous measurement must work for every state, and thus our analysis is 
sufficient to rule out these kinds of continuous measurements. (Further, as we will show 
below, measurements with Levy collapse operators tend to produce Gaussian states 
anyway.) Consequently, only the Gaussian (a = 2) case can correspond to a continuous 
measurement. 

3.2. Chained processes and exotic noise 

Despite the above result, it is possible to construct continuous measurements that 
contain the exotic statistics of any stable Levy process. To do so we combine the 
stable process with a Poisson process in the following way: we replace the time index 
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upon which the stable process L a (t) depends with the Poisson process V(t), to form 
the "chained" process S a (t) = L a (V(t)). (In the mathematics literature this chaining 
procedure is referred to as subordinating L a with V [ ].) Now examine the behavior 
of this new Levy process. Between jumps of the Poission process the time index for 
the Cauchy process does not change, and thus the increment dS a is zero. Upon an 
event (a jump) in the Poisson process, the time index of the stable process increases by 
unity, and thus generates a finite increment AL a = j Q dL a . The temporal scaling of 
the stable process L a no longer prevents the process S a from representing a continuous 
measurement, since the scaling as At — > is entirely determined by the Poisson process. 

When there are many jumps in a given time interval, the number of jumps fluctuates 
only a little from its mean value. Because of this, if we make the rate of jumps A large, 
and set a oc 1/A so that the effect of each jump on the measured system is small, then the 
Poisson-subordinated Levy measurement realizes what is effectively a quasi-continuous 
Levy measurement, with At = 1/A. Our simulations show that the measurement records 
of these measurements are qualitatively the same as those of the stable Levy processes 
(e.g., Fig 3). 

Measurements described by the Levy processes S a can be realized by using the 
imaging setup of Fig. I. The photodetection events are the jumps of the Poisson 
process, so we need each of these to correspond to a measurement with the result 
Ar = (X) + ^ALa for some 7. This is accomplished by choosing the aperture so that 
the square of the collapse operator is the density P(AL a ), with 7 set by scaling the 
width of P. In Fig. 2b we give the aperture profile T(£) for a Cauchy measurement. 
This prepares a state with a Cauchy probability density for position, and we display the 
Wigner function for this state in Fig. 6. 



4. Anomalous diffusion and continuous measurement 

4-1. Action of an arbitrary position measurement 

So far we have focused on the statistics of the measurement results. Now we turn to 
the question of generating anomalous diffusion in the dynamics of a quantum system. 
In particular, we address the question of what happens to the quantum state under the 
action of a Levy (or any other) spatial projection operator. We will assume a Gaussian 
initial state, with probability density 

p(x,*) = -^e-* 2 / 4 * 2 , (9) 

where for convenience we assume the state is centered at x — 0. As we are interested 
in the quasi-continuous — and thus weak- measurement — limit, we will assume that the 
width of the projection operator is much wider than the width of the state. Given 
a Hermitian collapse operator f2, we have the (unnormalized) measurement reduction 
p — > QpQ, or p(x,x) — > Q 2 (x)p(x,x) in the position representation. The operator 
Q 2 (x) is very broad, though not necessarily centered at x — 0. Given that it varies 
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slowly over the extent of the state, we can expand the operator to second order as 

Q 2 (x) = Q + fiiz + Vt 2 x 2 + 0(x 3 ), = Q e ax ~ bx2 + 0(x 3 ) (10) 

provided that a = Qi/flo and b = a 2 /2 + f^/lV After dropping normalization factors, 
the reduced state is 

n 2 (x)p(x) oc e ^-^\-* 2 l^ oc e -(*-M) a /*^ > (11) 

where 

2aa 2 a 

/i = TT46^' T= VTTib^ 2 ( J 

are the new mean and standard deviation, respectively, 

There are several points to note here. For a Gaussian collapse operator, for 
example, the b parameter is constant, since it is only the a parameter that controls 
the displacement of the Gaussian from x — 0. This reiterates the well-known 
result that a Gaussian measurement will reduce the uncertainty of a Gaussian state 
deterministically — i.e., independent of the measurement result. In particular, b = 
l/2a' 2 , where a' is the standard deviation of the measurement function. For any other 
form of the collapse operator, a particular measurement may increase or decrease the 
uncertainty. It is only on average that the uncertainty of the state decreases. Again, 
the b parameter, which represents a particular function of the local curvature O2 of the 
collapse function at (x), controls how the uncertainty changes due to the collapse. For 
example, for a long-tailed distribution such as the Cauchy, the uncertainty decreases 
if the measurement result is close to zero, but increases if the measurement result 
is very far from zero. Also since the a parameter controls the displacement of the 
measurement operator, it likewise controls the shift of the mean of the state in response 
to the measurement, as we see from fi oc a. However, the main point of this section 
is this: in this regime of broad measurement distributions, the action of the collapse 
preserves the Gaussian form of the quantum state, independent of the form of the 
measurement operator. In particular, even measurements of stable Levy distributions 
preserve the Gaussian form of the quantum state. 

4-2. Collapse to a Gaussian wave packet 

Now we extend the analysis above and prove a somewhat surprising result: 

Theorem: In the absence of any Hamiltonian evolution, a sequence of many repeated 
Levy measurements generates a Gaussian wave function — independent of the initial 
state — just as do Gaussian measurements. 

Proof: A sequence of N Levy measurements of X corresponds to multiplying 
the initial position density, P(X) = |^(x)| 2 , by a sequence of Levy densities. This 
means that the characteristic function of the resulting density is the convolution of 
the characteristic functions of N Levy densities. The characteristic functions of the 
symmetric a-stable densities, L a (x), are Xa(s) = exp[—a a \s\ a + ifis], where a is a width 
parameter and \i is the mean [9, 10]. For a > 1, \i is the mean of the Levy density. 
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Figure 6. Wigncr function for a state with a Cauchy probability density for position. 
Here X and P are scaled so that [X, P] = i, and the width of the Cauchy density is 
a = 0.3. Blue is positive, green negative. 



If we now think of Xa{s) as a probability density (for s), then convolving N of them 
corresponds to adding N random variables with density Xa- Since every nas a 

finite second moment, the central limit theorem tells us that for large N the result of 
these convolutions tends to a Gaussian. (Note that the usual central limit theorem does 
not handle complex-valued distributions, so we extend the theorem to cover this case 
in Appendix A.) If the final characteristic function is Gaussian, then the final density 
is also Gaussian. Thus, while a single Levy measurement on an initially flat density 
creates a Levy density, many repeated Levy measurements generate Gaussian densities, 
and thus Gaussian wave functions. □ 

4-3. Anomalous diffusion induced by back-action: a proposed experiment 

Once the wave-function of a single particle has been reduced to a Gaussian, even the 
exotic measurements discussed above induce only Gaussian noise in the dynamics of X. 
This is not difficult to show analytically for very strong (small 7) and very weak (large 
7) measurements, and the numerical simulation we have performed confirm it more 
generally. So how does one make measurements that do generate anomalous diffusion? 
We now show that one can do this using imaging, and in a manner that directly exploits 
the quantum nature of the measurement back-action. This distinguishes it from the Levy 
noise in Ref. [2] (as well as the momentum Levy flights in Refs. [17, 18, 19] that arise in 
subrecoil laser cooling), which can be treated semiclassically. After initially preparing 
the atom in a broad Gaussian state, we make a "square" measurement producing a wave 
function with a square profile; that is, a flat wave function with a sharp cut-off at each 
end. This is achieved by choosing the aperture to be a sine function. The quantum back- 
action of this measurement simultaneously generates a momentum probability density 
that is the square of the sine function. This has an infinite variance. We then wait for a 
quarter period of the harmonic evolution, and this transfers the momentum density to 
position. Repeating the "square-profile" measurement now produces an infinite-variance 
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time (arb. units) 

Figure 7. Eight trajectories of the mean atomic position, resulting from 
photodctection using an aperture with a square profile, and total efficiency of 50%. 

change in the mean position of the atom, and thus anomalous diffusion. Of course, in 
an experiment, each time we make a measurement, a number of photons will be lost 
before one is detected. We will refer to the average number of photons emitted for each 
one detected as the loss rate. We must check that the anomalous diffusion remains in 
the presence of this loss. 

We now simulate this experiment, assuming mirrors with a 130° arc. Using an 
aperture profile that gives a square-profile measurement with a width of 32 wavelengths 
results in a loss rate of 119. We note that the loss rate per measurement can be reduced 
by post-selecting those photodetections in which the photo was detected very soon after 
the laser pulse. For these detections there will be, on average, many fewer photons that 
were not detected. For example, in the present case post-selecting 1 in 118 measurements 
reduces the loss rate per post-selected measurement to unity. 

To show that our square-profile measurement induces anomalous diffusion in the 
mean position, (X), we must examine the statistics of the change to (X) caused by a 
single photodetection. Let us denote the change induced in (X) as A{X). To determine 
whether this change exhibits anomalous diffusion, we examine how the width of the 
probability density for A(X), which we will call o~i, compares to that for the sum of two 
independent samples of A(X), which we will call o~i. Anomalous diffusion is achieved 
when o"2 = 2 / Vi, with (3 > 0.5. To obtain accurate estimates of /3 we must simulate 
sequences with a large number of measurements. As a example, to determine (3 when 
the loss rate per measurement is 5, we simulated 24 sequences of 600 photodetections 
each. This is numerically intensive due to the long tails of the wave functions involved, 
and so we employ a parallel computer. For a loss rate of unity (50% efficiency) we find 
that p = 0.69 ± 0.16, and for a loss rate of 5 (17% efficiency) p = 0.60 ± 0.07. This 
shows us that while P does decrease with increasing loss, the signature of anomalous 
diffusion remains. This indicates that one can generate exotic noise for considerably 
higher loss rates, although confirming this numerically is prohibitive with our present 
computing resources. Examples of the evolution of the mean atomic position, for a loss 
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rate of 50%, are shown in Fig. 7. 
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A. Appendix: Modified Central Limit Theorem 

The usual version of the central limit theorem states that, given a random walk of 
independent, identically distributed steps, the resulting distribution converges to a 
Gaussian in the limit of many steps N : with a width scaling as y/N. The variance 
of each step must be finite for the central limit theorem to hold. Under certain other 
conditions, the one-step distributions need not be identical. To support our argument 
in Section 4.2, we will now modify the standard proof (see, e.g., Ref. [20]) of the central 
limit theorem, where the one-step distribution is the Levy characteristic function 

XaM = e~^ a+i ^ s , (A.l) 
where fi n is randomly chosen on each "step" from the corresponding Levy probability 
distribution P a {x) with zero mean. Obviously, this is not a sensible probability 
distribution. However, the central limit theorem examines the successive convolutions 
of the one-step distribution, which is well-defined, and ultimately what we are interested 
in. The characteristic function (Fourier transform) of the one-step distribution (A.l) is 

Xa,n(k) = 1 - ikfx n - fc2(ffw2 2 +/i " 2) + 0(k% (A.2) 

where we have introduced an overall scaling factor to normalize the one-step distribution, 
and the first two "cumulants" are given by 



fl n := / ds SXa,n( S ) 

J r (A.3) 
a n := / ds(s- /i„) 2 X Q ,n(s). 



(Recall that the terms in the expansion of the characteristic function give the successive 
moments of the original distribution.) In general, both of these quantities are nonzero, 
but the integrals are finite. Also, random variable with zero mean and finite 

variance. 

Now we will consider N "steps" taken in the random walk, but first it is convenient 
to rescale the width of the one-step distribution to make it narrower by a factor of 
\/N. This amounts to the replacements s — > sy/N and k — > k/\/N, so that the 
characteristic function of the characteristic function becomes 

X' a jk) = 1 - ^ - k2{ %^ l) + 0(N~n (A.4) 
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After A" steps in the "random walk," or the successive convolution of N of the one-step 
distributions, we have simply the product of the corresponding characteristic functions 
via the convolution theorem (again dropping overall factors): 

N r 



<*) = n 



n=l 



1 - 



ikpL r 



N 



2N 



0{N 



-3/2> 



The logarithm of the characteristic function of the convolution is 



A? 



n=l 

Expanding the logarithm gives 
logx« 

n= 

ik 



N 

£ 

71=1 L 



ikfl n k 2 (a 2 + ft 2 ) 



N 



2N 



+ 0(N- 



-3/21 



(A.5) 



(A.6) 



ikJjL n k 2 a 2 



0(N- 



-3/2 



N 



A- 5 



N 



(A.7) 



AT^''^ 2AT 

n=l n=l 



The usual central limit theorem applies to the first sum, which is O(yN), so the first 
term remains finite with unit probability as A" — > oo. In the same way, the stochastic 
part of the second sum vanishes as A" — > oo, and the second term becomes —k 2 (a 2 ) /2, 
where the angle brackets denote an expectation value with respect to fi n . Letting 
k — > —ik, we obtain the cumulant-generating function for the distribution after A^ 
steps. We see that in the limit of large N, the first two cumulants are finite and the 
rest vanish — this is the cumulant-generating function for a Gaussian distribution. 
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